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Bethe Ansatz solution of the stochastic process with 
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Abstract. The eigenfunctions and eigenvalues of the master-equation for zero range 
process on a ring are found exactly via the Bethe ansatz. The rates of particle exit 
from a site providing the Bethe ansatz applicability are shown to be expressed in terms 
of single parameter. Continuous variation of this parameter leads to the transition 
of driving diffusive system from positive to negative driving through purely diffusive 
point. The relation of the model with other Bethe ansatz solvable models is discussed. 
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The Bethe Ansatz pL^ is one of the most powerful tools to get exact results for 
the systems with many interacting degrees of freedom in low dimensions. The exact 
solutions of one- dimensional quantum spin chains and two-dimensional vertex models 
are classical examples of its application [2 . The last decade the Bethe ansatz has 
been shown to be useful to study one-dimensional stochastic processes The first 

and most explored example is the asymmetric simple exclusion process (ASEP), which 
serves as a testing ground for many concepts of the nonequilibrium statistical physics 0. 
Yet several other Bethe anstazt solvable models of nonequilibrium processes have been 
proposed such as multiparticle-hopping asymmetric diffusion model 0, generalizations 
of drop-push model [ZllHlin], and the asymmetric avalanche process (ASAP) [TUj . 

All these models has a common property. That is, a system evolves to the 
stationary state, where all the particle configurations occur with the same probability. 
This property can be easily understood from the structure of the Bethe eigenfunction. 
Indeed, the stationary state is given by the groundstate of evolution operator, which is 
the eigenfunction with zero eigenvalue and momentum. Such Bethe function does not 
depend on particle configuration at all and results in the equiprobable ensemble. Except 
for a few successful attempts to apply the Bethe ansatz to systems with nonuniform 
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stationary state, such as ASEP with blockage jTTj or defect particle ^21, there is still 
no much progress in this direction. 

In the other hand many interesting physical phenomena like condensation 
in nonequilibrium systems ^^1; boundary induced phase transitions ^1 or 
intermittent-continuous flow transition jTOl CSj become apparent from non-trivial form 
of stationary state, which change dramatically from one point of phase space to 
another. Typical example is the zero-range process (ZRP), served as a prototype of one- 
dimensional nonequilibrium system exhibiting the condensation transition [13,. While 
its stationary measure has been studied in detail ,17 , the full dynamical description is 
still absent. The aim of this Letter is to obtain the Bethe ansatz solution of ZRP. We 
obtain the eigenfunction and eigenvalues of the master equation of ZRP with special 
choice of rates, which is imposed by integrability requirement. We describe the phase 
diagram of the model and show that its particular cases are equivalent to ASEP, drop- 
push model and noninteracting diffusing particles. We also discuss the relation between 
ZRP and ASAP. 

Let us consider the system of p particles on a ring consisting of N sites. Every site 
can hold an integer number of particles. Every moment of time, one particle can leave 
any occupied site, hopping to the next site clockwise. The rate of hopping u^ni) depends 
only on the occupation number of the site of departure i. The stationary measure 
of such a process has been found to be a product measure [13], i.e. the probability 
of conflguration {rij} specifled by the occupation numbers {ni,n2, . . . yU^} up to the 
normalization factor is given by weight 

N 

W{M) = l[f{n,) (1) 

i=l 

where f{n) is given by 

/(^) = n ^ (2) 

if n > 0, and /(O) = 1. 

Another way to specify the conflguration of system is to use the set of coordinates 
of p particles {xi} = {xi, . . . , Xp}, rather than occupation numbers {ni, . . . , n^r}, the 
two ways being completely equivalent. Let us consider the probability Pt{xi, . . . , Xp) for 
p particles to occupy sites xi, . . . ,Xp at time t. It obeys the master equation deflned by 
the dynamics described. While the explicit form of the master equation is too bulky 
to be written down explicitly, particular examples will be given below to illustrate the 
solution. The main idea of its solution is to look for eigenfunction in the form of the 
Bethe ansatz weighted with the stationary weights of corresponding conflgurations (Q) 

Ptixu ...,xp) = W{{ni})Pt\xu ...,xp) (3) 

where Pt{xi, . . . , Xp) is the usual Bethe function deflned as follows 

{(Tl,...,(Tp} i = l 
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Here the summation is taken over all p\ permutations {cii, . . . , cr^} of the numbers 
{1, . . . ,p), and the coordinates of particles are assumed to be ordered in the increasing 
order xi < X2 < ■ ■ ■ < Xp. 

To explain the details we first consider two-particle case subsequently generalizing 
it to the case of arbitrary p. Without lost of generality we can define 

u{l) = l, u{2) = u>0. (5) 

Then we note that the function 

P°(Xi,X2) = e^* (Ai,2^[^4' + ^2,1^^4^) (6) 

satisfies the equation 

9iP°(xi,X2) = P°(xi - l,a;2) + Pi°(xi,X2 - 1) - 2P^{xuX2), (7) 

provided that 

X = zi + Z2-2. (8) 

When {x2 — Xi) > 2, the same equation describes the probability Pt{xi,X2) for ZRP. 
For xi = X2 — 1 = X the equation for ZRP 

dtPt{x, x + l) = Pt{x -l,x) + uPt{x, x) - 2Pt{x, x + 1) (9) 

can be obtained from[7j) if we define 

Ptix,x) = -Pt\x,x), (10) 
u 

which also follows from 0. It has been noted above that Pt{xi,X2) can be expressed 
trough P^{xi,X2) in the domain X2 > Xi. Therefore, the term P^{x,x — 1), which 
appears in the r.h.s. of (|7j) when X2 = Xi = x, does not carry any physical content. The 
freedom in its definition can be used to make the last equation for two-particle ZRP 

dtPt{x, x) = Pt{x- 1, x) - uPt{x, x) (11) 

consistent with ((Tj). If we redefine P^{x,x — 1) as follows 

P°(x, x-l) = {u- l)P°(x - 1, x) - (u - 2) P°(x, x), (12) 

((Tj) and (jlip become equivalent. The relation (jl2j) can be thought of as a constraint for 
the function P^{xi,X2)- After substitution of the ansatz © this constraint results in 
the relation between the amplitudes Ai^2 and ^2,1 

Al,2 ^ (2 - m) - (1 - u)z2 - Zi 

A2,l {2-u)-{l-u)zi- Z2 ^ ' 

Taking the cyclic boundary conditions, P°(xi, X2) = P°(a;2, Xi+N), into account reduces 
the problem to the system of two algebraic equations. The first one is the following 

^-N _ (2 - n) - (1 - u)z2 - zi 

{2-u)-{l-u)z,-Z2' ^ ' 

while the second can be obtained by the change Zi ^ Z2- 
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A generalization of these results to the case of arbitrary p leads us to the following 
constraint for -P°(xi, . . . , a;„) 

(nH-l)P°(...,(^-l),(^r\...)- 

n 

Y.Pt'{■■■,i^y-\i^-l),i^r\■■■)' (15) 

(«H-n)(...,«",...) = 

where (z)" denotes n successive arguments equal to i of the function P°(xi, . . . , x„). 
The Bethe ansatz to be applicable this relation should be reducible to the two particle 
constraint (|T2|l . This can be proved by induction. To this end, we assume that similar 
relations including u{k) are reducible for all k < n. Then starting from the relation 
fll5|) . which includes rate u{n), we apply (fT^ to every pair (z, z — 1) of arguments of 
the function P°(. . .) under the sum and require the result to be the similar relation 
including u{n — 1). We obtain the following recurrent formula for the rates 

u{n) = I ~ {I - u)u{n - I) (16) 

which can be solved in terms of g-numbers 

u{n) = H, (17) 



where 



1 _ 



with q = u — 1. 

Thus, the solution of the master equation for ZRP with special choice of the rates 
(I16II18P is given by the ansatz ()3I4|) with the eigenvalue 

X = j2z^-P, (19) 
i=l 

where the numbers Zi, {i = 1, . . . ,p), satisfy the system of algebraic equations. 



p 



(2 — m) — (1 — XL)Zj — 



The physics of the model can be extracted from the analysis of these equations. This 
study is beyond the goals of this Letter and will be published elsewhere. Nevertheless, we 
can get some qualitative information about the phase diagram of the model considering 
the behaviour of u[n) and particularly the limiting cases of parameter u reproducing 
the models known before. First, we note that for m = 2, i.e. g = 1, the g-numbers 
degenerate into simple numbers, \p\q=\ = n. Therefore, the rates are given by u{n) = n, 
which corresponds to the diffusion of noninteracting particles. The Bethe equtions in 
this case are reduced to 

zf = 1, (21) 

as is expected in noninteracting case. In the domain u > 2 the rates u{n) grow 
exponentially with n, which corresponds to the interaction between particles effectively 
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accelerating free diffusive motion, i.e. the higlier density of particles the faster is their 
mean velocity. In the limit u ^ oo the model is equivalent ton = 1 drop push model [7j, 
which is confirmed by the same form of Bethe equations |B]. In the domain 1 < u < 2 
the rates u{n) grow monotonously from u for n = 2 to 1/(2 — m) for n — > oo, resulting 
in the interaction between particles slowing down the particle flow. When m = 1, all the 
rates do not depend on the number of particles, i.e. u{n) = 1. This case can be mapped 
on the ASEP by insertion of one extra bond before every particle, 

{xi, . . . , XpjzRp ^ {xi + l,...,Xp + p}asep- (22) 

This mapping results in the totally asymmetric ASEP on the lattice consisting of N + p 
sites. At last, in the domain < m < 1 the rates u{n) also saturate to the constant 
l/{2 — u) with growth of n, though oscillating around this value. 

We should note, that both ASEP and drop-push model has been claimed to belong 
to Kardar-Parisi-Zhang (KPZ) [THj universality class|3llS!- Similarly, we expect that 
KPZ behaviour holds through all the driven-diffusive phases below and above the point 
u = 2 and this point itself corresponds to purely diffusive motion. The transition 
from one driven-diffusive (KPZ) phase to another with opposite driving through purely 
diffusive point has been studied in context of KPZ - Edwards- Wilkinson crossover in 
partially asymmetric ASEP 1^ . It would be interesting to compare the scaling 
behaviour with that appearing in our model. 

To show the formal relation of ZRP to another model solved by the Bethe ansatz, 
we first remind the formulation of the ASAP JUl ■ Consider particles at one dimensional 
periodic lattice. If all particles occupy different sites of the lattice, every particle cam 
jump to the right neighboring site with a Poissonian rate 1. When particle comes to 
already occupied site, the avalanche starts. It develops by discrete steps according to the 
following rules. If there are n > 1 particles at the same site, either all n particles move 
to the next right site at the next time step with probability /i„ or n — 1 particles move 
one step to the next right site with probability (1 — /i„), one particle remains unmoved. 
The avalanche is assumed to be infinitely fast comparing to the Poissonian time, i.e. 
it finishes before any new jumps of particles from sites other than active site. Now, 
let us look at an avalanche from the reference frame related to an active site. In this 
reference frame the particles from all sites other than active site move simultaneously to 
the left with probability 1. At every time step one particle from the active site occupied 
by n particles can also move to the left with probability 1 — fin- Obviously in this 
formulation the ASAP is nothing but the special case of discrete time ZRP. Indeed, the 
stationary measure of ASAP has been shown to be a product measure fSl- Moreover, 
the integrability imposes the same constraint on the toppling probabilities in ASAP as 
on the rates of exit from a site in ZRP studied in this Letter, which can be expressed 
through the relation u{n) = 1 — /i„ for u < 1. It is interesting that the solution of 
ASAP was based on usual form of the Bethe ansatz, resulting in the uniform stationary 
measure. It, however, allowed the consideration of the stable configurations of particles 
only, while avalanches played the role of nonlocal interaction and led to the master 
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equation containing infinite series. 

Despite the formal similarity, there is a big difference between ASAP and ZRP 
proposed above. Note, that in the ASAP the dynamics allows at most one active site to 
exist at any time moment. This situation corresponds to the inhomogeneous version of 
ZRP with one attracting site. Evans ^Hj has shown that in this case the condensation 
transition always takes place. So it does in the case of ASAP, corresponding to 
the transition from the intermittent to continuous flow phase. In the other hand, 
ZRP studied in present Letter is spatially homogenous, while the rates u{n) saturate 
exponentially fast to the constant with growth of n when u < 2. Hence, they do 
not satisfy Evans condensation criteria, which requires u{n) to decay slower than 2/n. 
However, the large fluctuations for small u are expected, leading to phenomenon like 
sharp crossover. It is of interest to see how such fluctuations affect the relaxation 
dynamics. This is a matter for further investigation of the Bethe ansatz solution. 

To summarize, we used the Bethe ansatz weighted with the stationary weights 
of particle configurations to solve the stochastic process with nontrivial stationary 
state. We found the eigenvalues and eigenfunctions of the master equation for zero 
range process in continuous time. The system of Bethe equations was obtained. We 
showed that the Bethe ansatz to be applicable, the rates of exit of particle from a site 
should be expressed in terms of q-numbers depending on single parameter. For special 
values of this parameter the model corresponds to the drop-push model, ASEP, and 
noninteracting particles. The variation of this parameter was argued to lead to the 
transition similar to KPZ - Edwards- Wilkinson crossover in the partially asymmetric 
exclusion process. The relation of the model with asymmetric avalanche process was 
discussed. 
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